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Abstract 

Extending earlier work of R. Donaghey and P. J. Cameron, we investigate some canon- 
ical "eigen-sequences" associated with transformations of integer sequences. Several known 
sequences appear in a new setting: for instance the sequences (such as 1, 3, 11, 49, 257, 1531, 
. . .) studied by T. Tsuzuku, H. O. Foulkes and A. Kerber in connection with multiply transi- 
tive groups are eigen-sequences for the binomial transform. Many interesting new sequences 
also arise, such as 1, 1, 2, 26, 152, 1144, . . ., which shifts one place left when transformed by 
the Stirling numbers of the second kind, and whose exponential generating function satisfies 
A'(x) = A(e x - 1) + 1. 
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1. Transformations of sequences 

The purpose of this note is to present certain integer sequences that are associated in a 
canonical way with various sequence transformations. This is an extension of earlier work 



of R. Donaghey [13|, who was concerned with the transformation denoted by BINOMIAL in 
the present paper, and of P. J. Cameron ||, who mainly considered the transforms that we 
call EULER, INVERT and WEIGH. Cameron says that the last of these transformations was 
suggested to him by Professor Seidel. We hope that Professor Seidel will enjoy the sequences 
presented here. 

A number of papers have appeared in recent years in the combinatorial literature dealing 



with generalizations of the Catalan and Bell numbers (e.g. [14|, [16|, [49|, [p6| ). The present 
paper provides another way to view certain of these sequences, as well as suggesting still further 
generalizations . 

We shall not give any proofs of the identification of certain of our sequences with those in 
the literature — in most cases these are easily obtained from known properties of the sequences. 
The majority of the sequences discussed here are not in [54], although they will appear in [ ]55j , 

A typical sequence will be written a = [ao, a\, . . .], with ordinary generating function (o.g.f.) 

oo oo 

A(x) = &nX n and exponential generating function (e.g.f.) A{x) = J2 a n x n /n\. The 

n=0 n=0 

logarithmic derivative A'(x)/A(x) is denoted by £{x). Unless indicated otherwise, "sequence" 
means "sequence of real numbers". The following list contains the transformations that we 
shall consider. Here a is the original sequence and b the transformed sequence. In each case b 
is integral if a is. 

The initial transformations map a sequence a = [oq, ai, . . .] to b = [bo, bi, . . .]. 
BINOMIAL 

bn = £(fja k , a n = ±(-l)^(fjb k , (1) 



with two equivalent forms 



B(x) = -±-A , B{x) = e*A(x) . (2) 



References: J46|, p. 4], |20|, p 192], |2l|, p. 6], |13|], ]16| (in the latter reference the inverse of 
this transformation appears under the name "Euler transformation" ) . If a n enumerates a class 
of structures on n points in which every element is involved, b n enumerates the same type of 
structure when only a subset of the points needs to be involved. 

This transformation also arises when analyzing sequences by means of difference tables. 
The familiar difference table associated with a sequence [co,ci,...] is its difference table of 
depth 1. By using the leading diagonal of this table as the top row of a new table, we obtain 
the difference table of depth 2, and so on. 

For example the difference tables of depths 1 and 2 for [1, 3, 9, 27, . . .] are 

1 3 9 27 1 2 4 8 

2 6 18 12 4 

4 12 12 
8 1 

To recover the sequence from the leading diagonal of the difference table of depth r, we apply 
the binomial transform r times. In the example we go from [1, 1, 1,1,.. .] to [1, 2, 4, 8, . . .] to 
[1,3,9,27,...]. 



STIRLING 



or equivalently 



b n = S(n, k)a k , a n = ^ s(n, k)b k 
k=0 k=0 



B(x) = A(e x - 1) , (3) 



where S(n,k) and s(n,k) are respectively Stirling numbers of the second and first kinds (the 
numbers of the second kind being the more fundamental) fl46| , pp. 90,202] p0| , p. 296], pl| , 
p. 7]. If a n is the number of objects in some class with points labeled 1,2, ... ,n (with all labels 
distinct, i.e. ordinary labeled structures), then b n is the number of objects with points labeled 
1,2, ... ,n (with repetitions allowed). 

The following transformations are related to the operation of convolution. They do not in 
general have inverses. 
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CONV 



EXP-CONV 



F-CONV 



b n = ^2 akdn-k, B(x) = A(x) 2 . 

k=0 



bn = E ( k ) a kan-k, &{x)=A(xf 



k=0 



k=0 

where F(r, s) is one of LCM, GCD, AND, OR, XOR (in the last three F is applied to the 
binary representation of its arguments). The F-CONV transforms were suggested to us by M. 
Le Brun (3^], and will be applied only to sequences of nonnegative integers. 
The remaining transformations map a = [01, 02, . . .] to b = [b±, 62, . . .]. 



MOBIUS 



b n = ^2n (j^j ad, a n = ^2b d 



d\n d\n 

where jx is the Mobius function. Two equivalent forms are 



00 00 7 

, n s , n s 

n=l n=l 
n=l n=l 



the latter being a Lambert series [27, p. 257]. 
WEIGH 

00 00 

i + Y J b n x n = n(i+x n r . 

n=l n=l 

6 n is the number of ways of getting a weight of n grams if we have aj weights of % grams, i > 1 



1 44, p. 174]. Cameron [0] gives other interpretations. 



EULER 

00 00 -. 

1 + y b nX n = tt — . 

n=l n=l " > 

If a n enumerates a class of connected structures on n unlabeled nodes, b n enumerates the same 
structures when connectedness is not required. b n is also the number of ways of partitioning 
the integer n, given that there are possible types of parts of size k, for k = 1,2,.... 
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Alternatively, if a n is the number of generators of degree n of a graded polynomial algebra, 
then b n is the dimension of the n-th homogeneous component, i.e. the number of linearly 
independent monomials of degree n. WEIGH has a similar interpretation in terms of graded 



exterior algebras. References Q, [51] and || give other interpretations. Calculations are 

oo n— 1 

facilitated by setting the left side equal to exp E c n x n so that c n — nb n — E c^6 n _^., 

d\n 

PARTITION. For 1 < a\ < a% < . . ., form c = \c%, C2, ■ ■ ■] by discarding all duplicate entries 
from a, so that c is strictly increasing; then 

oo oo 
n=l n=l 



b n is the number of ways of partitioning n into parts that can have sizes ai, 02, 



INVERT 



x 



n=l 1 - E °n^ n 

n=l 

b n is the number of ordered arrangements of postage stamps of total value n that can be formed 



if we have types of stamps of value i, i > 1 |44, p. 174]. || gives several other interpretations. 



00 00 
REVERT. If y = E a n x n with a x = 1, define b by x = £ (-l)" +1 &„y" (so 61 = 1) @, 

n=l n=l 

p. 149], 0, p. 16], @]. 
EXP 

1 + B(x) = exp^(x), A(x) = log(l + B(x)) . 
This transformation relates the number connected structures on n labeled nodes to the total 



number of structures 0] , [ 51 1 . 
2. Eigen-sequences 

Let T denote any of the above transformations. We are interested in eigen-sequences a 
such that T o a = a, or (since in most cases that equation has only trivial solutions) 

T o a = Loo, T 2 o a = L o a, T o a = L 2 o a , 

T o a = N o a, T o a = M o a , 
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etc., where the auxiliary operators R, L, M, N are defined by 

Ro[a ,ai,a 2 ,...} = [1, a , ai, . . .] , 

Lo [a ,a 1 ,a 2 ,. ■ ■] = [a 1} a 2 , a 3 , . . .] , 

N o [a ,ai,a 2 , . . .} = [00,-01,-02,...], 

Mo [00,01,02,...] = [oi,2oi,2a 2 , • • •] . 

Note that the e.g.f.'s for Roa and Loa are respectively l+Jo* «4.(z)d:c and .4.' (x). The equations 
T o a = L o a and T o a = M o a were studied by Cameron H] for the transformations EULER, 
INVERT and WEIGH. 

Suppose T is linear, in the sense that it maps a to b with 

n 



k=0 



Let -D denote the array [L> nifc : n = 0, . . . ; k = 0, . . .], and set D r = [D%[]. Then a sequence 
satisfying T r o a = L s o a has the recurrence 

n 

a n+s = D ni a k . ( 5 ) 

A:=0 

and, provided D n)7l = 1 for all n, sequences satisfying Toa = NoaovToa = Moa have 
recurrences 

-y n — 1 

On = -- J! D n,kO,k , (6) 

or 

n-l 

On = X! ^",fc a fc ' (7) 
fc=0 

respectively. 

Tables I and II show the most interesting of these eigen-sequences. The symbol in the last 
column of Table I indicates the sense in which the sequence is unique: 

a: unique sequence fixed by the operator; furthermore any sequence converges to this one 
under repeated application of the operator 

j3: unique sequence beginning 1, . . . that is fixed by the operator 

7: unique sequence [0,1,02,03,...] of nonnegative integers that is mapped to [0,02,03,...] 
under XOR-CONV 

5: every sequence converges to one of these five (see below for further details) 



e: lexicographically earliest integer sequence [1 = a\ < < ■ ■ •] fixed by REVERT. 

We now comment on some of the individual sequences, especially those that have already 
appeared in the literature. 



BINOMIAL. SI is the familiar sequence of Bell numbers, with A(x) = e eX , C(x) 



|3J, |pq| , IM , |10|, 210]. a n is the number of ways of partitioning a set of n objects (Motzkin's 
"sets of sets" pip . This sequence has property a, as indicated in the table. In particular, 
BINOMIAL o a = L o a. In fact b = L o a = [1, 2, 5, 15, 52, . . .] is the unique sequence b that 
begins 1, . . . and has the property that the leading diagonal of its difference table is Rob. Al- 
though this property of the Bell numbers must be well-known, we have not seen it mentioned 
before. a n is also the number of orbits on n-tuples of an n-fold transitive group. Sequences 
S2, S3, S4 (and further sequences of the same type) appear in Kerber [||] (see also Tsuzuku 
1 59], Foulkes ]l8|]) in connection with the characterization of multiply transitive permutation 
groups by cycle structure. These sequences have C{x) = e 2x , e 3x , e 4x , etc., a pleasant property 



not recorded in [29]. Furthermore b = L o S2 is the unique sequence b that begins 1, . . . and 
has the property that the leading diagonal of its difference table of depth 2 is R o b. Similar 
properties hold for the other sequences. Sequence S6 has ^4(x) = 1/(2 — e x ) (from Eq. @), 
and gives the number of unlabeled planar rooted trees with n + 1 terminal nodes all at the 
same height, in which at every other height there is a node with at least two successors Q, 
1 39]. E.g. d2 = 3 refers to 

A\ 





Equivalently, a n is the number of preferential arrangements or total preorders of n things 



(Motzkin's "lists of sets"), i.e. the number of linear orderings of n things using < and = [22|, 
|4l[| . The orderings corresponding to the above trees are 1<2, 2 < 1, 1 = 2. Sequence S49 has 
A(x) = 2/(1 + e x ) and is related to the Bernoulli numbers B n ]l|, 810], the tangent numbers 



T n |3l[] , and the H n numbers of Terrill and Terrill |57| by 

2 '->n+l 



n + 1 



(2" +i - l)B n+1 (n>0), 

• T„. ^(n+l)/2 



= = -qsgy* („>i) . 

v ' n 2( n+1 )/ 2 ~ 

Donaghey [O] has studied a different kind of self-inverse sequence associated with this trans- 
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form, namely sequences satisfying a n = (~ 1)^0^™ ka ki where d = 25 is a constant. 

k=0 

STIRLING. Only one of these sequences seems to have been studied before. For sequence 

n-1 

S10 we have A(e x — 1) = 2A(x) — 1, a n = J2 S(n,k)ak (n > 1). This sequences gives the 

fc=o 



number of rooted (but not planar) trees with n + 1 labeled nodes all at the same height [34] 
E.g. a2 = 4 refers to 






1 o 2 



Equivalently, a n counts a certain class of preferential arrangements on unordered pairs of 
elements of an (n + l)-set, called "ultradissimilarity relations" |52|| , 
The e.g.f.'s for S7 and S50 satisfy 

A'(x) - A(e x - 1) = 1 , 
A{x) +A{e x - 1) = 2 

respectively. Recurrences for all these sequences follow from (|)-(|). 
The identity 

([0, p. 209]) shows that 

STIRLING o R = R o BINOMIAL o STIRLING , 

(with a little care about the zeroth term), from which it follows that STIRLING maps the 
all-ones sequence to the Bell numbers (SI). 

CONV. Sll is the ubiquitous sequence of Catalan numbers, one realization of which is the 
number of unlabeled planar rooted trees in which each node has < 2 successors [|l(], p. 53], 



|4|, ||. Sequences Sll, S12, S13, etc. satisfy 

CONV r o a = L o a , 
hence xA(x) r = A(x) — 1, and are given by 

1 / rn\ 
n (r — l)n + 1 \ n J 



These sequences have a long history pOl Section 7.5, Example 5], p0|. They also appear in 



the context of enumerating clusters of polygons in 24]. 
The sequences satisfying 

CONVoa = L s oa, (s = 2,3,...) 

appear to be new. The first of them, S14, has 

, , . 1 - y/l - 4x 2 - 4x 3 
A{X) = 2^ ' 

Except for S15, the factorial numbers, all the sequences emerging from the other convolutional 
transforms are new. 

MOBIUS. From here on the subscripts of the sequences start at 1. 
S23, with a\ = 1, 

a n +i = ^2 a d, n > 1 , 

d\n 



is the number of planted achiral trees [19], |pq |. 
S26, with a\ = 1 and 

a n=Yl a <i' ( n > 2 ) . 



d\ n 
d < n 



is the number of ordered factorizations of n. Its Dirichlet generating function T>(s) = a n/ nS 

n=l 

is 2/(1 + Q(s)~ l ). a n is also the number of "perfect partitions" of n — 1, that is, partitions of 
Ti — l which contain precisely one partition of every smaller number (i.e. a set of weights such 



that every weight of m grams, m < n — 1, can be realized in exactly one way) |p7| |, [45, p. 123] 



Tg, p. 126]. 

The other sequences appear to be new. Sequence S51 has Dirichlet generating function 
V(2 -C(s)- 1 ). 

WEIGH. S27 gives the number of unlabeled rooted trees with no symmetries (|26|], see in 
particular Eq. (1*)). That this sequence shifts left under WEIGH was already pointed out by 
Cameron ||. 

EULER. S30 counts unlabeled rooted trees, and satisfies 

OO OO -| 

j2 a n +ix n = n — 



n=0 



li 1 (1 - X n Y 

n=l v ' 



|45| , p. 127], p5| . S33 counts series-reduced planted trees, and Mo 533 counts series-parallel 
networks ^?J, p. 142], |36[| . The connections of these three sequences with the EULER 
transform have been investigated by Cameron ||, 



PARTITION. Every sequence converges to one of these five. S34 is fixed, while S35 swaps 
with S36, and S37 with S38. There are no other sequences of finite order under this trans- 
formation. In particular, S34 is the unique sequence [01,02,...] with the property that the 
number of partitions of n into parts of sizes {01, 02, . . .} is equal to a n . We find it remarkable 
that, in spite of all the papers that have been written about partitions, these five sequences 
seem never to have been published before. 

INVERT. Cameron ||] has already observed that the Catalan numbers (Sll) shift left under 
INVERT, and that S39, the solution to Schroder's second problem satisfies 
INVERT o a = M o a. Like the Catalan numbers, the latter sequence occurs in many guises, 



sometimes being called super-Catalan numbers |p3 |, fig , p. 168], |32j] , |LCj, p. 57], fl49| , [p5 | 



|56], [23]. a n is among other things the number of rooted planar trees with n + 1 nodes. The 
o.g.f. satisfies 2A(x) 2 — (1 + x)A(x) + x = 0, which implies the inversion property mentioned 
above, as well as the fact that this sequence also satisfies INVERT 2 o a = L o a. 

S41 is the sequence of Motzkin numbers, again a widely-studied sequence |4C], R, [16|, 
|l|, p|. The o.g.f. satisfies A(x) 2 = (1 + x)(A(x) - x). 



Sequences S42, S43 and similar sequences are related to the Motzkin numbers, and have 



been studied in pC|], [56], pg|], in part because they enumerate secondary structures of nucleic 
acids. 

REVERT. S44 is new, and no combinatorial setting is presently known for it. 



EXP. S45 is the sequence of Euler numbers, with e.g.f. (1 + sin x)/ cos a; [43, p. 110], |12| , 
p. 262], H p. 259]. 

S48 is the solution to Schroder's fourth problem, which is also connected with series-parallel 
networks g|, p. 197], |K| p. 224], |§. The e.g.f. satisfies exp^4(x) = 2A{x) + l-x 0, p. 224, 
Eq. (aa)]. 

S52 arises in an unexpected context. If we multiply a n by 2™" 1 we obtain an integer 



sequence [1, —1, 1, 1, —13,47, 73, . . .], sequence 2093 in [54]. This occurs in numerical analysis 
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in the computation of Airey's converging factor for asymptotic series, both for the Weber 
parabolic cylinder function and for the exponential integral ||, [38], fnfl , [42|. The simple 
equation for the e.g.f. for S52 which follows from the eigen-sequence property, 



exp.A(x) = l + 2x- A{x) , 

does not seem to be given in any of these references. Once again the other sequences are new. 
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Table 1(a) Integer eigensequences 



# Sequence Operator Property 



51 


[1 


1 


2 


5,15,52,203,877,4140,...] 


R o BINOMIAL 


a 


52 


[1. 


1, 


3. 


11,49,257,1539, 10299,...] 


R o BINOMIAL 2 


a 


53 


[1. 


1. 


4. 


19,109,742,5815,51193,...] 


R o BINOMIAL 3 


a 


54 


1> 


1 


5. 


29 201 1657 15821 1 


R o BINOMIAL 4 


Qj 


55 


[1, 


1, 


1. 


2,4,9,23,65,199,654,...] 


R 2 o BINOMIAL 


a 


56 


[1, 


1, 


3. 


13,75,541,4683,47293,...] 


M- 1 o BINOMIAL 


P 


S7 


[1, 


1. 


2. 


6,26,152,1144,10742,...] 


R o STIRLING 


a 


S8 


[1, 


1, 


3. 


14,97, 934,11814,188650,... 


R o STIRLING 2 


a 


S9 


[1, 


1, 


1. 


2,5,16,66,343,2167,...] 


i? 2 o STIRLING 


a 


S10 


[1, 


1, 


4. 


32,436,9012,262760,...] 


M~ l o STIRLING 


P 


51 1 

kJ J. J. 


[1 


1 


2. 


R 1 4 42 1 32 429 1 430 1 


R o CONV 


r\ 
LX 


512 


[1. 


1. 


4. 


22,140,969,7084,53820,...] 


R o CONV 2 


a 


513 


[1. 


1, 


8. 


92,1240, 18278,285384,...] 


R o CONV 3 


a 


514 


[1. 


1. 


1. 


2,3,6,11,22,44,90,187,...] 


R 2 o CONV 


a 


515 


[1. 


1, 


2. 


6,24,120,720,5040,...] 


o EXP-CONV 


a 


516 


[1. 


1, 


4. 


28,280,3640,58240,...] 


R o EXP-CONV 2 


a 


517 


[1. 


1. 


1. 


2,4,10,30,100,380, 1600,.. 


i? 2 o EXP-CONV 


a 


518 


[1, 


1, 


2. 


5,14,40,128,369,1214,...] 


R o LCM-CONV 


a 


519 


[1, 


1, 


2. 


3,4,6,6,11,10,18,16,...] 


R o GCD-CONV 


a 


520 


[1. 


1, 


2. 


1,2,4,0,5,2,4,0,10,0,...] 


o AND-CONV 


a 


521 


[1, 


1. 


2. 


7,20,58,174,519,1550,...] 


R o OR-CONV 


a 


522 


[o, 


1. 


2. 


4, 14,38,118,338, 1006,...] 


i2 o XOR-CONV 


7 


523 


[1. 


1. 


2. 


3,5,6,10,11,16, 19,26,...] 


R o MOBIUS^ 1 


a 


524 


[1. 


1, 


3. 


5,10,12,24,26,43,52,78,../ 


R o mobius~ 2 


a 


525 


[1, 


1. 


1, 


2,2,4,3,7,4,11,6,15,7,...] 


R 2 o MOBIUS" 1 


a 


526 


[1, 


1, 


1. 


2,1,3,1,4,2,3,1,8,1,3,3,... 


] M- 1 o MOBIUS" 1 


P 
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Table 1(b) Integer eigensequences (cont.) 

# Sequence Operator Property 



527 


[1, 


1 


1 


2,3,6,12,25,52,113,247,...] 


R o WEIGH 


a 


528 


[1, 


1 


1 


3,6,16,43,120,339,985,...] 


R o WEIGH 2 


a 


529 


[1, 


1. 


1, 


1,2,2,4,6,10,17,29,51,...] 


R 2 o WEIGH 


a 


530 


[1, 


1, 


2 


4,9,20,48,115,286,719,...] 


R o EULER 


a 


531 


[1, 


1. 


3 


8,25,77, 258,871,3049,...] 


R o EULER 2 


a 


532 


[1, 


1. 


1. 


2,3,6,10,20,36,72, 137, ...1 


R 2 o EULER 


a 


533 


[1, 


1. 


2 


5,12,33,90,261,766,2312,...] 


M" 1 o EULER 


/? 


534 


[1, 


2. 


2 


4,5,7,9,12,16,20,25,32,...] 


PARTITION 




535 


[1, 


2. 


2 


4,4,6,7,11,12,16,18,25,...] 


PARTITION 2 




536 


[1, 


2. 


2 


4,4,7,8, 12, 13, 18,21,29, ...] 


PARTITION 2 




537 


[1, 


2. 


3 


4,6,9, 11, 15, 19,25,31,41, ...] 


PARTITION 


6 


538 


[1, 


2. 


3 


5,6,10,12,17,22,29,36,48,...] 


PARTITION 2 


6 


511 


[1, 


1. 


2 


5,14,42,132,429,1430,...] 


R o INVERT 


a 


539 


[1, 


1. 


3 


11,45,197,903,4279,20793,... 


] R o INVERT 2 


a 


540 


[1, 


1. 


4 


19,100,562,3304,20071,...] 


R o INVERT 3 


a 


541 


[1, 


1. 


1, 


2,4,9,21,51,127, 323,835,...] 


R 2 o INVERT 


a 


542 


[1, 


1. 


1, 


1,2,4,8,17,37,82,185,423,...] 


R 3 o INVERT 


a 


543 


[1, 


1. 


1, 


1,1,2,4,8,16,33,69,146,...] 


i? 4 o INVERT 


a 


539 


[1, 


1. 


3 


11,45,197,903,4279,20793,... 


] M' l o INVERT 




544 


[1, 


2. 


4 


7,10,12,18,40,44,45,...] 


REVERT 


€ 


515 


[1, 


1, 


2 


6,24,120,720,5040,...] 


iioEXP 


a 


545 


[1, 


1, 


1. 


2,5,16,61,272,1385,...] 


R 2 o EXP 


a 


546 


[1, 


1. 


1, 


1,2,5,15,53,213,961,...] 


R 3 o EXP 


a 


547 


[1, 


1. 


3 


14,89,716,6967,79524,...] 


R o EXP 2 


a 


548 


[1, 


1. 


4 


26,236,2752,39208,...] 


M _1 o EXP 


P 
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Table II Fractional eigensequences 



# 



Sequence 



Operator 



Property 



549 
550 
551 
552 



1 -1 - ■ 



2 ' u ' 8 ' ' 



31 n 691 



4 ' 



111 

2' 4' 2' 



19 



39 2623 365 
16' 32 ' 4 ' 



' 2 ' 2 ' 4' 2 ' u > 2 ' 8 ' 4 ' u > 2 ' 8>2> u '"- 



111 

"2' 4' 8' 



2' 8' 4' 

13 47 73 
16 ' 32 ' 64 ' 



2' 8' 2' 

2447 16811 
128 ' 256 ' 



N- 1 o BINOMIAL p 

N- 1 o STIRLING 

N- 1 o MOBIUS /3 

TV" 1 o EXP p 
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